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Entanglement of superconducting qubits via microwave fields: classical and quantum 

regimes 

Jian Li, K. ChalapatQ and G. S. ParaoaniQ 

Low Temperature Laboratory, Helsinki University of Technology, P.O. Box 5100, FIN-02015 TKK, Finland 

We study analytically and numerically the problem of two qubits with fixed coupling irradiated 
with quantum or classical fields. In the classical case, we derive an effective Hamiltonian and 
describe its entangling properties. We identify a coupling/decoupling switching protocol and we 
construct composite pulse sequences leading to a CNOT gate. In the quantum case, we show that 
qubit-qubit-photon multiparticle entanglement and maximally entangled two-qubit states can be 
obtained by driving the system at very low powers (one quanta of excitation). Our results can be 
applied to a variety of systems of two superconducting qubits coupled to resonators. 

PACS numbers: 03.67.Lx,85.25.Cp,74.50.+r 



I. INTRODUCTION 

In recent years, there has been a consistent experi- 
mental progress in the quantum-coherent manipulation 
of superconducting circuits based on the Josephson ef- 
fect. Singlc-qubit operations have been demonstrated by 
now in charge flux Q, charge-flux 0], and phase [3| 
qubits. Some of these experiments [|[ made it clear that 
in order to achieve long decoherence times the qubits 
have to be operated at the so-called optimal point, where 
the first order noise induced by fluctuations of external 
control parameters (gate voltages, magnetic fluxes) can- 
cels. The notable exception from this rule is the phase 
qubit, which by construction does not have an optimal 
point; this is compensated by engineering a peculiar bias 
circuit 

The first-generation experiments with fixed coupling 
0, 0] did manage to achieve two-qubit gates but the 
qubits were not operated at the optimal points, with a 
corresponding loss in fidelity. Later it was noticed that 
flux qubits can be operated at the optimal point if the 
coupling is realized through a dc-SQUID modulated at 
the sum and difference of the qubits' resonance frequen- 
cies 01- If the dc-SQUID is replaced by a third flux qubit, 
largely detuned from the frequencies of the two qubits, 
the coupling can be realized through the quantum induc- 
tance of the additional qubit Q . 

With the advent of circuit QED architectures 0, [l(| , 
several qubits can be placed in the gap between the signal 
line and the ground of a coplanar waveguide resonator. In 
the dispersive qubit-resonator coupling regime, the qubit- 
qubit coupling are realized by virtually exciting the res- 
onator. Furthermore, due to the structured, cleaner elec- 
tromagnetic environment around the qubits, relatively 
long decoherence time can be achieved. 

Finally, since any additional coupling elements tend to 
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introduce extra decoherence in the system, one can place 
a further restriction and ask the question: is it possible 
to devise schemes in which the qubits are operated at 
the optimal points and there are no active additional ele- 
ments for coupling? The first such proposal was inspired 
from NMR: the so-called FLICFORQ protocol [lj, with 
coupling realized through the dressed states of each qubit 
under on-resonancc microwave fields. The protocol was 
extended to more general detuned driving fields [l2| , but 
although theoretically sound, it has not been realized yet 
experimentally. In the FLICFORQ protocol, the fixed 
coupling is much smaller than the qubit-qubit detun- 
ing; even under the enhancement the effective coupling 
strength is up to only one eighth of the bare coupling 
strength [13, [H| • A small coupling strength is desirable 
for singlc-qubit operations, but limit the two-qubit gate 
speed. 

Therefore, one very interesting question emerges: is 
it possible to develop a protocol which retains the ad- 
vantages of FLICFORQ and has also the capability of 
performing fast two-qubit gates? In this work we will 
analyze such a protocol, in which the two qubits are rel- 
atively strong coupled in the absence of driving fields; 
moreover, we show explicitly how to construct microwave 
pulses to switch the coupling off. 

The outline of this paper is as follows. In Sec. [II] wc 
study a specific circuit consisting of two qubits coupled 
to a single cavity mode (a coplanar waveguide resonator). 
A general fixed qubit-qubit coupling Hamiltonian is de- 
rived. We then analyze the case in which the qubits are 
manipulated using classical fields: in Sees. IHII and HVl by 
deriving the effective Hamiltonian of the system in the 
rotating reference frame, we show that for large values 
of the driving field amplitudes, the coupling is switch- 
able. We also find numerically the on/off coupling ratio 
of the switchable coupling, as well as the validity of our 
rotating wave approximation (RWA). The implementa- 
tions of a single-qubit gate and a CNOT gate are demon- 
strated in Sec. [V] We show that both gates can be re- 
alized with high speed and high fidelity. In the quantum 
case, in which the qubits interact with a single quanta of 
radiation, we propose a quantum nondemolition (QND) 
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method to entangle the qubits based on measuring the 
presence of an excitation in the resonant cavity to entan- 
gle the qubits in Sec. IVI| and finally we conclude our 
work in Sec. IVIIl 



II. MODEL HAMILTONIANS 

Although many of the results derived in this paper 
are rather general (once the Hamiltonian is put into any 
of the forms used in this paper), it is useful to start 
by analyzing a concrete superconducting quantum cir- 
cuit, which will serve as our workhorse. Consider a 
system of two coupled charge qubits, irradiated with 
monochromatic off-resonance microwave fields. The cir- 
cuit is shown in Fig. [TJ Two split single Cooper pair box 
(SCB) qubits are capacitively coupled to the center con- 
ductor of a coplanar waveguide resonator. Each box has 
its own read-out circuitry (such as a large current-biased 
Josephson junction, as in the case of charge- flux qubits). 




m± 



FIG. 1: Schematic circuit of our model system. 

To ensure first-order insensitivity to charge fluctua- 
tions, we voltage-bias the gates of the SCBs at the opti- 
mal points 



CglVdcl — Cg2Vi c 2 — 



(1) 



When two microwave fields with the same angular fre- 
quency ujd and different phases 0x,2 are applied, the gate 
voltages of the SCBs have AC components 



= V^ w i,2(*) cos(uj d t + 0i, 2 )- 



(2) 



Written in the eigenbasis of each SCB | f) = (|0) + 
|l))A/2 and | |) = (|0) - |l»A/2, |0) and |1) denoting 
two lowest charge states, the Hamiltonian has the follow- 
ing form (see Appendix 



hscb — 

3=1,2 



E 



E C jWj(t) COs(u! d t + <t)j)(Tj 



(3) 

where Ejj and E@j indicate the standard Josephson and 
charging energies, and Wj(t) = —C g jV^ W j(t)/2e. 



A large current-biased Josephson junction (CBJJ) is 
also coupled to the end of the resonator through a ca- 
pacitor; the bias current is such that, during many of 
the qubit operations described in this paper, only two 
bound states (of energy difference uj\q close to that of 
the resonator) are relevant. This junction will serve as a 
detector of the state of the resonator, as we will see in 
Sec. IVI1 The Hamiltonian of the CBJJ assumes the form 
of a two-level system, 



H 



CBJJ — 



UJlO 
<7 



(4) 



where u>iq denotes the transition frequency between the 
two lowest levels in the tilted cosine potential. 

To simplify the equations, from now on our units will 
be such that H = 1. 

Depending on the length of the stripline, the qubit- 
qubit and qubit-CBJJ couplings must be treated in dif- 
ferent ways. For a relatively long resonator, the couplings 
are mediated by excitations in the resonator. As derived 
in Appendix [Al the total Hamiltonian in this case is 



Hscb + Hcbjj + uj r (^a + 1/2) 

+i gj [a* + 2iVj(t) cos(w d i + <f>j)] (a* - a) 

3=1,2 

+ti{a+ - (7 3 ")(a t - a), 



(5) 



with oj r the resonance frequency of the resonator, (a) 
the photon creation (annihilation) operator, gj (n) the 
qubit-resonator (CBJJ-resonator) coupling strength, and 
a t ( a ;T) the CBJJ raising (lowering) operator. By per- 
forming a RWA to neglect fast oscillating terms, we have 



H 



totl 



SCB 



+ H cb jj +u r {(Ja+ 1/2) 



3=1,2 

— n(a^a + 'a'), 



(6) 



where o\-~ = | J,)(f | and a J = \ | denote qubit rais- 
ing and lowering operators, respectively. Here we have 
neglected the terms describing interactions between the 
qubit driving fields and the resonator (in the second line 
of Eq. {5)). This is justified by the fact that Rabi fre- 
quencies satisfy = EcjWj(t) <C Ejj, and the SCBs 
are working in the charge regime Ecj 3> Ejj, therefore 
Wj (t) should be much smaller than unity. 
a. Long resonator limit 

In the long resonator limit (the resonator length of the 
same order as the wavelength corresponding to an exci- 
tation with energy of the order of qubit energy) , a fixed 
linear qubit-qubit coupling can be achieved in the disper- 



sive regime, gj <C Aj 



- Ejj . By performing a gener- 



alized Schricffcr- Wolff transformation (see Appendix [C]) 



exp 



3=1,2 3 



aaf) 



3 



\- g j n j {i)cos{u} d t + <j) ] ) f 



on the Hamiltonian ([6]), we obtain 



H 



totl 



E 



^-cr| + % COs(uJ d t + <\>j)o X ; 



5ig 2 (Ai + A 2 ) 



(7) 



2A!A 2 

Here we assume that the CB J J is biased far off resonance 
with the resonator such that it is effectively decoupled 
from it. 

b. Short resonator limit 

For a very short resonator (the resonator length much 
smaller than the wavelength corresponding to an excita- 
tion with energy of the order of qubit energy), the center 
conductor can be considered as a small metallic island. 
The two qubits and the CBJJ are capacitively coupled 
to this island. In Appendix [B] we show that the total 
Hamiltonian has the form 



At, 



SCB 

E 



H 



CBJJ 



'l u 2 



E 

i=i,2 



a j a 3 



U> 2 (£) cos(uj d t + q 
+wi(t) cos(uJdt 4 
ljWj(j)cos(uJdt 



3 • 



(8) 



The Hamiltonians © ^ and dTTJ) are similar, there- 
fore we will use a generic time-dependent Hamiltonian 



H(t) 



E 

J=l,2 



■j-a] + Qj cos(uj d t + (f>j)a^ 



(12) 

to replace all of them, uij' and u> xx denote the Larmor fre- 
quency of qubit-j and the qubit-qubit coupling strength, 
respectively. As discussed in Sec. HI we are interested in 



the regime u> a 



1.2- 



III. SWITCHABLE COUPLING MECHANISM 

Now, we start to derive an effective Hamiltonian with 
tunable coupling terms. In order to eliminate the ex- 
plicit time dependence, we bring the Hamiltonian Q12p 
into the rotating reference frame by transforming it with 
the operator 



S\(t) = cxp 



iuo d t 



K + ^2 



(13) 



By performing a RWA to neglect oscillating terms 
with frequency 2uj d , we get a time- independent effective 
Hamiltonian 



ffeff = H x 



E 

j=l,2 



(14) 



with £a 2 /4 and "fj/2 the qubit-qubit and qubit-CBJJ 
coupling strengths, respecti vely . Due to small coupling 
E12 <C Eci.2, the crosstalk [12j terms in the second and 
third lines of Eq. ([5]) are negligible. We can also drop 
the interactions between the qubit driving fields and the 
CBJJ, described by the last line in Eq. (|5|), because Wj is 
small. Therefore the total Hamiltonian is approximately 

Hscb + Hcbjj 



H 



tot2 



E W„x„x \ ~~ 7j 



I (9) 



If the coupling capacitances C m \ t2 
coupling is dominating 



> a 



m3; 



the direct 



E 



E 

i=i,2 
, E12 



+ flj COs(LU d t + 4>j)<7j 



l u 2- 



(10) 



For C m i )2 <C C m 3, the two qubits can be coupled through 
virtual excitation of the CBJJ. By considering dispersive 
coupling 7j <C A'j = t^io — Ejj , the qubit-CBJJ couplings 
are eliminated by another Schrieffer- Wolff transformation 
similar to the one in long resonator case, and we obtain 



H 



tot2 



E 

i=i,2 



E 



Jj 



+ flj cos(u; d t + (t)j)<7j 



7172 ( A^ + A' 2 
4A' : A 2 



■(0i"o- 2 



■)■ (11) 



where 



Hj = + Y^ cos ^ a 3 ~ sill( ?W) 



is the Hamiltonian for qubit-j, and 



H 12 = —(aM + afaV) 



(15) 



(16) 



describes the interaction between the two qubits. 6j = 
u> d — ujL denotes the detuning between the driving fre- 
quency and the corresponding qubit Larmor frequency. 

By diagonalizing the single-qubit Hamiltonians given 
by Eq. (fl~5j) . we obtain the eigenenergies and the corre- 
sponding eigenstates of each single qubit: 



E 



U) 



-55j/2, 



)U) =cos ^ e ^72i 



T) 



sin -j-e-^/ 2 1 |), (17) 



sin ^ e ^/2| t)+cos ±e-^l 2 \ 1), (18) 



where 



^/5 2 + O 2 , sin0 7 - = cos0 7 = (19) 



We then project the time-independent effective Hamil- 
tonian (TTJ]) onto the new basis of product states | + +) , 
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■+)> 



), and get 



H 'cS 



r\ Z 1 o z 



.XX 



X 

XX 



3S(^)| sin^iCT^ 1 ' + cos6>icr^ 

n jo-W sin 6» 2 ct< 2) + cos 6 2 a { x 2) 

>}, (20) 



sin^crW 



(T 



(2 



with the Pauli matrices crS 1 ^ and (s — x, y, z) in the 
new basis, and the phase difference <fi = <pi — <p 2 . 

Next wc transform H' cff into a new rotating frame with 

(21) 



S 2 (t)=expj-^ 



.(2) 



By assuming that Wx,2 are much larger than the coupling 
strength uj xx , we may perform a second RWA to neglect 
oscillating terms with frequencies toi. 2 and u>i + uj 2 . The 
resulting Hamiltonian is 



H'Ja 



— cos 0(1 + cos 6>i cos 6 2 



) jcos(<5wi) 



^(1)^(2)4.^(1) (2) 
w x x 1 w y " y 



sin </>(cos 6»i + cos 6 2 ) {sin(5arf) cri 1 ^ 2 ' + <4 1)(T y 2) 



- sin(<5wt) 
cos((So;i) 



^(1)^(2) _ a (D CT (2) 
w y x w x w y 



a {i) a m _ a {i) a (2) 

u y u x w x u y 



■ cos <p sin t/i sin #20z c 



(1)^(2) 



(22) 



with Suj = uji — uj 2 . This effective Hamiltonian has oscil- 
lating terms. In the rest of this section, we discuss how to 
switch the effective coupling off in two limits with respect 
to the oscillation frequency Suj. 



A. High 8uj 

The oscillating terms in (f2"2")l can be neglected if |<5o;| ^> 
us xx . By using the definition of sin^- in (TiU)l . the remain- 
ing effective Hamiltonian is expressed as 



■ cos< 



Of) 



(23) 



A maximum coupling strength of about lo xx /2 can be 
achieved when the driving frequency is in-resonance with 



the qubit Larmor frequencies u>d 



Ji 2 , and large driv- 
are applied in the 



ing amplitudes fii 2 ^ \Si, 2 \ ~ 
meantime. 

To turn the coupling off, we can switch off either J7i or 
il 2 . Nevertheless, the two conditions must be satisfied: 
uji- 2 > uj xx and \8uj\ > uj xx . Without loss of generality, 
we assume Vl 2 = 0. The condition uii i2 ^> lo xx is fulfilled 
when the driving frequency is largely detuned from the 
qubit frequencies, Si ~ S 2 ^> oj xx (assuming ujd > wf 2 )- 
The second condition is fulfilled by driving the first qubit 
with a rather large amplitude fix 3> \/2Siuj xx . 

We notice that this decoupling mechanism can also be 
derived in the case of quantized driving fields: let us 
start with the usual Jaynes-Cummings form [l3| of the 



Hamiltonian 



H = uidO^ a — 



J=l,2 



3 _z 



>l u 2 



g(cr^a + a x a^), 
(24) 

and assume a dispersive coupling g <C 5\. Then we per- 
form a Schrieffer- Wolff transformation U = cxp [g(cr^ a — 
<Jia))/5i\ to eliminate the direct qubit-field coupling to 
leading order, and obtain 



U^HU 



LUdCL a 



UJ 2 



9_ 

Si 



. .XX —X —X 

■ OJ a 1 a 2 ■ 



(25) 



Due to the driving field applied on it, the Larmor fre- 
quency of qubit- 1 is ac-Stark shifted by the quantity 
2g 2 (a { a) / 5\. If the ac-Stark shift is much larger than 



the qubit-qubit coupling strength ui xx 
arc effectively decoupled. 



the two qubits 



B. Low Suj 

A more general situation is when Suj is not so large, 
therefore the oscillating terms can not be neglected. 
Since we are considering a relatively small qubit-qubit 
detuning, it is not possible to eliminate 1 + cos^i cos 6 2l 
cos 0i + cos 6 2 , sin 61 sin# 2 , and satisfy uj\^ 2 3> ui xx at the 
same time. However, there is still a way to switch the 
coupling off. As shown in Eq. (|2"2"|) . terms in the first 
and the third lines have a common factor cos d>. If the 
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second line can be removed by setting 8j and ilj to obtain 
Si 8 2 



COS Vi + COS V-2 



0. 



(26) 



the rest of the Hamiltonian will be switchable by means 
of </>. One solution of (|26[) is 5i = #2 = 0, which is realized 
only when the qubits are on resonance. For off-resonance 
qubits, Eq. (|2"6"|) leads to 



(27) 

By defining A = 
TUl) we obtain the 



An extreme case is when 8u> = 0. 
uji —UJ2, 81 = 82 — A, and using Eq. 1 
expression of the resonance condition 



nl = nj + 28 2 a - a 



and the coupling coefficients 



a a 25f + Ql-5 2 A 

f + COS &i COS tf 2 = 



COS Vi + COS 6>2 



8l + nl 

28 2 - A 



sm 6*1 sin 6^2 



^2 + 28 2 A - A 2 



(28) 

(29) 
(30) 

(31) 



To eliminate cos6*i + cos6* 2 , we set the driving fre- 



quency LO d 



K 



y )/2. The resonance condition 



(|28p. as well as the condition 
fti = fi 2 > ui xx . Since |A| « uj X: 
nian has a rather simple form 



Si, 2 3> becomes 
: , the effective Hamilto- 



(32) 



Does this switching scheme work with quantized fields 
as well? Here we briefly notice that in this approach it 
is essential to operate with states of the electromagnetic 
field having well-defined phases; therefore this switching 
scheme can be implemented either with classical fields or 
with coherent states [l3|, UH; Fock (number) states, even if 
they could be prepared experimentally, have fluctuating 
phases, therefore cannot be used. 



IV. ENTANGLING PROPERTIES 

In order to characterize the effectiveness of the switch- 
able coupling schemes derived in previous section, wc 
need to determine the on/off ratio of the coupling with 
the original Hamiltonian (fT2"|) . In this section, we will use 
the concurrence [iH | to study the entanglement between 
the two qubits. The concurrence of a pure two-qubit state 
\ip) is defined as 

C^) = Mcj v ®cj y \r)\, (33) 
where \ifi*} is the complex conjugate of \ip). For a general 
two-qubit state \ip) = c uu \ H) + c ud \ ||) + c du \ |f) + 
Cdd\ II), it is 



C(V>) = 2\c uu c d d - c ud c d u\ < 1. 



(34) 



We define the on/off ratio of coupling as the ratio of 
the maximum concurrence obtained when the coupling is 
switched on to the maximum concurrence obtained when 
the coupling is effectively off. 

In the lab frame, the concurrence can be calculated by 
numerically solving the Schrodinger equation with the 
time-dependent Hamiltonian (fT2]) [l7l ]. In the rotating 
frame, analytical calculations with the effective Hamil- 
tonian (|14p can be done under the circumstance that 
81 and 82 arc small. Here wc only focus on the zero 
8uj case discussed in Sec. IIIIBI For the sake of sim- 
plicity, we also consider that the driving fields satisfy 
ill = ^2 = ^ 3> uj xx , and specify our initial state to the 
ground state | ft). 

We start by deriving the analytical expression of the 
concurrence. When 4>i ~ (f>2 = 0, the coupling is switched 
on. The effective Hamiltonian (TTJJ is approximately 



O 



(35) 



It has eigenvalues Ai 
± ^J{lo xx Y 



0, A 2 = — uj xx , A 3j 4 = 
± ^ {10** )* + 4il 2 ]/2, and the corresponding (un- 
normalized) eigenvectors [— 1,0, 0,1] T , [0, — 1,1, 0] T , 
[1, \ 3 ^/£l, A 3j 4/r2, 1] T , respectively. By expanding | ff) 
in terms of the eigenvectors, we get the time-dependent 
state vector 



\m) 



l(e-^ + e -^ + 2) + _^!_ 



'g — iXit _ e -i\3t^ 



TT) 



-iX^t ^ — iX^f) 



2(A 3 - A 4 ) 
1 



(ITD + IIT)) 



4^ 



iX^t 1 — 1X4 



e~ lM t-2) 



4(A 3 - A 4 ) 



II) 



(36) 



r 



Then the concurrence is approximately 



oscillating between and 1 with a period of 2ir/ui x 



C(t) « \e iu * - l|/2 = y/[l - cos(cj**t)]/2, (37) 
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According to (|32|) . if we set fa = n/2 and fa — 0, the 
coupling should be effectively switched off, whereas the 
effective Hamiltonian (fl4l) in this case reads 



The maximum concurrence is about (2w xx ) 2 / XI 2 . 



(a* 2 



-(of of + 



(38) 



By diagonalizing H c ff and expanding | f"f) with the eigen- 
vectors again, we arrive at an analytical expression of the 
time-dependent concurrence as 



C(t) 



(u xx ) 2 
2il 2 



|3 + cos(2fii) -4cos(Oi)|. (39) 



Thus, the on/off coupling ratio is approximately 
il 2 1 (2ui xx ) 2 . To perform the numerical calculations, we 
set Q = 10uj xx , cu d = 200w xa; , = uj d + A/2, and 
ll>2 = — A/2 (considering A <C u)d)- For a charge 
qubit (Ec/Ej » 1), the probability of leakage to non- 
computational states is negligible for such a ratio of 
Q/u>i 2 (sec Appendix |D|) . 



(a) i)A = Oco xx ii) A = 2 <o xx (b) i)A = 0co xx ii) A = 2 co' 




FIG. 2: (Color online) Time evolution of the concurrences for the input state jV'in) = | TT)- ( a ) — 4>2 = 0; (b) 4>i = n/2, 
<j>2 = 0. The solid (red) and dashed (blue) lines indicate the concurrences calculated with the time-dependent Hamiltonian (|12p 
and the effective Hamiltonian (|14[) , respectively. We have set Q = Wcj xx , cod = 20f2, and uii t 2 = ujd ± A/2. 



In Fig. [2l we show the concurrences calculated with 
different values of A. The numerical results indicate that 
the switchable coupling scheme works well as long as A 
is not too big compared to uo xx . An on/off coupling ratio 
of about 20 can be obtained in a time duration longer 
than several 7r-pulse widths. Once A becomes large, 
the on/off ratio is reduced. Furthermore, differences be- 
tween the concurrences (dashed lines) calculated with the 
effective Hamiltonian (fT4)) and the exact concurrences 
(solid lines) evaluated with the time-dependent Hamil- 
tonian (fT2"]) start to appear. These differences imply that 
the counter-rotating terms neglected in the RWA leading 
to (IT4l) cause effects larger than the Bloch-Siegert shift 



If 



In order to understand how the counter-rotating terms 



Hr 



U! xx atcrt 



j=l,2 



.#3 , 



D -2iui d t 



(40) 



3=1,2 



behave, we extend the derivation of Bloch-Siegert shift 
for single two-level atom (see e.g. Chapter 7 of Ref. [3), 
to our two-qubit system. To shorten the following dis- 
cussions, we only consider = fa = 0. 

The time-evolution operator generated by (fT2")) can be 
expressed as 



U(t) = S*i(t)e^ rfTt ^exp -if drH a (r) 

l Jo 



(41) 



where Si(t) is dehncd in Eq. (fHT]) . i/ cff is given by (fl"4")l . 
'T exp[- • •] denotes the time-ordered exponential integra- 
tion, and 



JH eB T 



H rr e 



-4H e ffT 



(42) 



Because of the qubit-qubit coupling terms in H e fi , the 
transformation (|4"2"|) is hard to perform. Due to the fact 
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that only when Si 2 = TA/2 3> oj xx the deviations be- 
come significant, we may ignore the coupling terms, and 
carry out (42|) in the basis of {| ++), | +-), | -+), | )} 



H ci (t) 



cx Pi ~ 



x cxp 



H'. 



IO\0 



(i) 



(2) 



. (43) 



where lo\ = u>2 = \J (A/2) 2 + Q 2 , and H' CI denotes the 
counter-rotating terms transformed into this basis. 




FIG. 3: (Color online) Time evolution of concurrences un- 
der the same values of parameters as those in Fig. [2ja iv). 
Again, the solid (red) and dashed (blue) lines indicate the 
exact concurrence and the concurrence evaluated by the ef- 
fective Hamiltonian with RWA. The dotted (green) line is the 
concurrence calculated with counter-rotating terms compen- 
sation. 

We have evaluated the time-ordered exponential inte- 
gration to second order. By transforming it back to the 
computational basis {| ||),| |J.),| |t)>l we have- 

obtained 



T 



exp 



1 1 dri? cr (r) 
Jo 

I- if drH c ,(T)- f dr 2 f 2 dTiH a (T 2 )H CI ( n ) 
Jo Jo Jo 



exp 



LO d UJ\ 



A 



\Pl ~ G 2 J 



(44) 



This equation shows that besides the Bloch-Sicgert 
shift in each qubit Larmor frequency, the counter- 
rotating terms also give rise to Rabi frequency shift pro- 
portional to A. As illustrated in Fig. [31 the concurrence 
recalculated by substituting (|44f in (|4Tj) (dotted line) fits 
the exact concurrence (solid line) much better than the 
concurrence calculated barely with H c r (dashed line). 
V. IMPLEMENTATION OF QUANTUM GATES 
BY COMPOSITE PULSES 



We are now in the position to address the problem of 
implementing quantum gates. Simple rectangular driving 
pulses have been proved to be not satisfactory for precise 
manipulations [21j . It is better to adopt other pulse tech- 
niques such as shaped pulses and composite pulses. 
In order to achieve arbitrary quantum gates with both 
high accuracy and high speed, we need to find a numer- 
ical approach to get the suitable pulse parameters (am- 
plitudes, frequencies, phases and pulse durations). The 
idea is to find a certain functional, and to obtain the 
optimal pulse parameters by maximizing or minimizing 
this functional. In a previous work [23| , we found a set of 
parameters to get a locally equivalent gate of CNOT by 
minimizing the functional |Gi(t)| 2 + IG2 (£) — 1| 2 , where 
G\ and Go, are two locally invariant functionals defined by 
Makhlin [2J|. This method has a drawback that, when 
a locally equivalent gate is found, determining the two 
singlc-qubit gates which transform it into CNOT could 
be a difficult task. Therefore in this section, we minimize 
a different functional (2f| which can lead us directly to 
the target gate we want. 



We work in the rotating frame to reduce the computa- 
tional complexity. Since in the preceding section we have 
found that the RWA broke down when A was large, here 
we assume A = 2ui xx , and set the detunings Si and 82 
to be —A/2 and A/2, respectively. The pulse duration 
of performing a gate, t p , is equally divided into m small 
intervals. We consider cj>i and 02 as "global" parameters 
which arc unchanged in the whole pulse duration. In 
order to avoid sharp edges as in rectangular pulses, we 
assume that in each time interval dt the "local" pulse pa- 
rameters Qi 2 vary linearly with time, and at t = and 
t = t p , r2i ; 2 = 0. The unitary operator generated by the 



effective Hamiltonian (|14| right after the pulses can be 
well approximated as 



U, 



eS 



exp < —idtH c fi 



(m-l) ^( m ) 



(m-1) 



■ exp 



-idtH, 



eft 



(0) (1) 



)(0) 



(45) 



where > 0) stand for the values of 0^2 at the end of the kth interval, and denote fh,2 at t = 



(0) 
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The optimal 0i.2, fij] and t p for a target gate U arc 
achieved by searching for the global minimum of the error 
functional 

e = yjTr[{U-U cS y(U-U cS )]. (46) 



A. Simulated annealing 

Since normally e has many local minima, to avoid be- 
ing trapped in a local minimum, we use the simulated 
annealing (SA) method [26L [27| to minimize e. It em- 
ploys a random search of pulse parameters which al- 
lows changes not only decreasing the "energy" e, but also 
increasing it. The probability to accept a change is 
given by P = exp(-<Je/T), where 6e = e a f tcr the change - 
^before the change, and T is a parameter acting as the sys- 
tem "temperature". One can easily find that, if a change 
of pulse parameters decreases e, P is always larger than 
unity, which means we always accept this change; if a 
change increases e, we still have chance to accept this 
change. Our algorithm is summarized as follows. 

(i) Define the starting temperature T s = — 0.01/ ln(0.8) 
and the halting temperature = — 10~ 10 / ln(0.8) for the 
annealing process. 

(ii) Define the boundaries of pulse parameters € 
[0, 10u> xx ] and 0i,2 € [-7r/2,7r/2]. 

(iii) Initialize the pulse parameters fij 2 = ^> ujXX an d 
0i,2 = 0. Calculate the initial value of e with these pulse 
parameters. Initialize the temperature T = T s . 

(iv) Repeat the following steps if T > 

a. randomly generate sequences of pulse parameters 
until all the parameters in the sequence are inside the 
boundaries; 

b. evaluate e with the last sequence; 

c. keep the resulting e and the pulse sequence with a 
probability P = cxp(— Se/T); 

d. after every 1000 successful evaluations of e, lower 
the temperature by 1%. 

We demonstrate two examples here: a 7r-rotation of 
the first qubit around its A-axis 



U Xl = e 



h = 



and a CNOT gate 



C^cnot — e 



— 27T/4 



-i 

-i 

-i 

-i 



10 
10 
1 
10 



(47) 



(48) 



where the factor exp(— in/A) is used to make C/cnot G 
SU(4), since Tr(H cff ) = 0. 

Because of the off-resonance shaped pulses and the fi- 
nite on/off coupling ratio, a pulse duration much longer 



than 7r/max(fi) = 0.lir/ui xx was expected for X\ gate. 
By running our SA program numerous times, we have 
found it very difficult to minimize e to a satisfactory 
value if t p < 0Att/oj xx . In Fig. 0Ja), a possible control 
sequence of {f2i(f), S^2 (*)} is presented. The pulse dura- 
tion t p = 0Att/ijj xx has been divided into 10 intervals, 
and the optimized phases <f>i = 0, 4>2 — — 0.57T. 

The minimum pulse duration for the CNOT gate we 
have found so far is t p = 1.2ir/u> xx . Interestingly, it is 
approximately equal to the summation of the interaction 
cost [H C ff (CNOT) = 0.5tt /uj xx for our effective Hamil- 
tonian, and durations for two singlc-qubit gates, which 
indicates that the SA program constructs the CNOT gate 
in a way close to a Cartan decomposition (29[. A possi- 
ble sequence of {f2i(t), ^(i)} is shown in Fig. BJb). The 
corresponding phases cf>i = — 0.2n and 4>2 = 0.077T. 
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FIG. 4: Possible sequences of fii(£) and Q2 (t) (a) for the 
single-qubit gate Xi, (b) for the CNOT gate. 
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B. Obtaining maximally entangled two-qubit states 




FIG. 5: (Color online) Trajectories of the reduced density 
matrices in the rotating reference frame, projected onto the 
Bloch spheres. The blue and the red dots denote the motions 
of qubit-1 and qubit-2, respectively. 

To get a visual picture of how these optimized driv- 
ing pulses work, we perform a simulation of producing a 
maximally entangled two-qubit state by the CNOT gate. 
Suppose that initially the two qubits were in the ground 
state | tt) and then a 7r/2-rotation around Y"-axis of 
qubit-1 has been applied; so the new input state is then 

W = M 8 | t ), (49) 

as indicated with arrows in Fig. 

We then send the qubits through a CNOT gate, re- 
alized using the pulse sequences presented in Fig. [4jb); 
then the output state of the two qubits will be a Bell 
state [3(ij 

l0n> = ^(ITT> + IU». (50) 

The total density matrix for this system of two coupled 
qubits is obtained by numerically solving the Schrodingcr 
equation with the effective Hamiltonian (fLE)) ; the reduced 
density matrix of each single qubit is obtained by par- 
tially tracing out qubit-1 or qubit-2. Fig. Oshows the mo- 
tions of the reduced density matrices in the Bloch sphere 
picture. The reduced density matrices of each qubit end 
up in the centers of each of the spheres (corresponding to 
a zero Bloch vector), indicating that the two qubits are 
in a maximally entangled state [30j . 

For the total density matrix of the output state in ro- 
tating frame p ou t, we find a state fidelity [22j 

Wtt^/w) = yfiMp^M > 0.999. (51) 



1 




0.98 



-0.1 -0.05 0.05 0.1 

uncertainty in co xx 

FIG. 6: The fidelity loss due to uncertainty in measuring lj xx . 

C. Losses of gate fidelity 

Now, we move back to the lab frame to evaluate the 
gate fidelity [H[ 

HU) = P^W%j (52) 

with the time-dependent Hamiltonian ([12"]) . We set 
0Jd = 200oj xx and 2 = ojd — 5\,2, as in Sec. |lVl By 
numerically solving the Schrodingcr equation, we obtain 
F{X X ) = 0.9993 and .F(CNOT) = 0.9982. Such high fi- 
delities can only be obtained in rather ideal cases, since 
the calculations have only taken gate errors due to the 
counter- rotating terms into account, and these errors are 
very small by choosing small qubit-qubit detuning. In 
practice, uncertainties in system parameters, external 
noises, etc. will cause extra fidelity losses. 

The optimization method we have used relies on the 
well defined qubit parameters, such as the Larmor fre- 
quencies u>i2 and coupling strength ui xx . In experiments 
on Josephson qubits, although the qubit parameters may 
be tunable due to specific designs, we still assume they 
arc static, and the detailed knowledge of them are ob- 
tained by measurements. So, uncertainties in measure- 
ments of qubit Larmor frequencies and coupling strength 
will give rise to gate errors. By performing numerical sim- 
ulations, we have found that, compared with JT(CNOT), 
F(Xi) is less sensitive to the uncertainties in qubit pa- 
rameters. As shown in Fig. [6l there is nearly no loss of 
^(Xi) even when the measured value of ui xx is 90% (or 
110%) of the exact ui xx , whereas in order to keep high 
.F(CNOT) (say > 0.99), the uncertainty in uj xx should 
not be more than 7%. The data in Fig. [7] indicate that 
the uncertainties in measuring u>i 2 should be controlled 
within 0.05%, which is easily achievable with present-day 
electronics. If we consider typical charge qubits and dis- 
persive coupling through a resonator, io\ 2 ~ 27r x 5 GHz 
and uj xx w 2tt x 20 MHz (see (3), the allowed unccr- 
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taintics arc in MHz range. 



F(CNOT) 




-2-1 1 2 

uncertainty in x 1 o 



FIG. 7: (Color online) Loss of ^(CNOT) due to the uncer- 
tainties in measuring 0J\2- ^(Xi) has similar features to 
.F(CNOT), but less sensitive to the uncertainties (data not 
shown) . 

Nevertheless, the uncertainties in the driving pulses 
are not troublesome to the optimized gates. As shown in 
Sec. IV Al the required precision of pulse parameters are 
already low, only 0.1lu xx for Qi ; 2 and 0.017T for <fii,2- To 
account for experimental imperfections (jitter and phase 
noise of the external electronics used to create the pulses) , 
we also simulate the errors in pulses as normally dis- 
tributed random numbers with relative error a and null 
average, and find that if a is smaller than 1%, the loss of 
fidelity is negligible. 

We now examine the effect of the electromagnetic de- 
grees of freedom which inevitably couple to each of the 
two qubits, producing decoherence. We consider here the 
worst-case scenario [201 ] . in which each qubit is coupled 
to a different environment, modeled by Lindblad super- 
operators 

A J [p] = ^-(a>|-p), (53) 

and 

AiW = —(2aTpa+-afaJp-pafaj), (54) 

describing longitudinal and transversal noise with decay 
rates and respectively. We evolve the system 
numerically under the simultaneous action of decoher- 
ence and of the pulse sequences corresponding to single 
gates and two-qubit gates. To simplify the presentation 
we take the decoherence rates of the two qubits equal, 
= y, T$ = Tfi, and we show the resulting fidelity 
loss due to decoherence in Fig. [S] 




T (or r^) / (B xx 

FIG. 8: Loss of two-qubit (CNOT) fidelity ^(CNOT) and 
single-qubit (X\) fidelity J-(Xi) due to decoherence. 



D. Summary 

Compared to the FLICFORQ protocol developed in 
[Til ] . our protocol do improve the speed of two-qubit gate, 
due to the fact that in our schemes the qubit-qubit cou- 
pling strength is reduced only by a factor of 2, and it 
still holds the advantages of FLICFORQ: no need to dc 
bias away from the optimal points and no need for active 
tunable coupling. For a relatively weak coupling medi- 
ated through cavity, lu xx ps 2tt x 20 MHz, a CNOT gate 
can still be performed in about 30 ns. With two qubits 
directly coupled by a capacitor, lu xx can be much larger, 
and it is only constrained by oj xx <C r2 12 "C 2 . 

We have only demonstrated quantum gates with the 
switchable scheme for low 8u> (see Sec. IIII B[) . however, 
the numerical optimization method used in this section 
should be applicable also to the high 5ui case (in Sec. 
IIII A[) . with one more pulse parameter needed to be 
optimized. Since the high Sui scheme can be realized by 
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singlc-qubit driving, there might be possible applications 
other than manipulating qubits, like bringing a qubit on 
and off resonance with a naturally-formed two level sys- 
tem (TLS) located in the Josephson junction barrier, as 
briefly discussed in Appendix [El 



resonator, has the form 



H = 



-igi 
-ig 2 



(55) 



VI. QUANTUM EFFECTS: QND 
ENTANGLEMENT 



In the preceding sections, we considered that the cou- 
pling strengths between the qubits and the resonator are 
much smaller than the detunings between them. Two- 
qubit entanglement can be realized by using the CNOT 
gate. In this section, we move to the non-dispersive 
regime, and propose a protocol to create maximally en- 
tangled two-qubit states using a method [H| inspired 
from cavity QED systems. 

We still bias the dc gate voltages of the qubits at charge 
degeneracy points. After exciting one qubit, we switch off 
the driving fields. The Hamiltonian of this qubit-qubit- 
resonator system is then 



H 



E 

.7=1,2 



E 



erf + u) r a) a + i 



E 

.7 = 1,2 



")■ 



It has eigenvalues 



Ai,2 = ±ysi + fla, A 3 = 0, 
with corresponding eigenvectors 



(56) 



|Ai, 2 ) 
|Aa) 



1 



(<?i|l> + <? 2 |2» 



<? 2 |l)+<?i|2» 



(57) 



For symmetric couplings g\ = g 2 = g, the eigenstate | A3) 
is a direct product of the resonator vacuum state and a 
maximally entangled two-qubit state. 

An arbitrary initial state \ipin) = ai|l) + a 2 |2) + aa|3) 



|ai| 2 + |a 2 | 2 + \a>3 | 2 = 1) can be rewritten as 



ijax +a 2 ) _a3_ 
2 V2 



By considering that the two qubits and the resonator 
have the same energy Eji = Ej% = & r , the Hamiltonian 
projected onto the basis states |1) = | J,f) g) |0) p , |2) = 
I T4> ® |0) p and |3) = I TT> ® |1) P , with |0) p and |l) p the 

zero-photon and one-photon Fock states of the field in and the time evolution of it 



|Ai> 



ai - a 2 



V2 



A3) 



i(ai + a 2 ) a.3 

2 + 7f 



|A 2 ), 



(58) 



-iXit 



i(ai + a 2 ) _a3_ 

2 + y/2. 



\Xi)+e 



-i\ 2 t 



i(ai+a 2 ) a 3 



|A 2 ) - e 7= — I A 3 ) 



V2 



a\ — a,2 1 ( at + a 2 ia 3 
2 2 



p -iV2gt , I ( fl l+ a 2 ^3 \ iV 2gt 



Q 2 - Q l , 1 / fl l +°2 _ ^3_A ^-JV^g* ,l( fl l + Q 2 , ^3 \ p iV2 9 t 

2 + 2>v 2 V2J 2 V 2 + v/2' 



|2) 



i(ai + a 2 ) 03 
2 



D -iV2gt 



2^2 

c 1 (t)|l)+c a (t)|2)+c 3 (t)|3). 



i(oi + q 2 ) 03 
2^ 2 



|3) 



(59) 



The global entanglement among the two qubits and 
one photon can be quantified by the Q-measure [3~ij | 

Q(t) - I [|ci(i)c 2 (t)| 2 + | Cl (i)c 3 (£)| 2 + \c 2 (t)c 3 (t)\ 2 ] , 

(60) 

which can be analytically calculated with Eq. ([59| . In 
Fig. numerical calculations of the Q value are illus- 
trated. 

Now, we consider the consequences of the fact that 



the resonator is lossy: we look at the case in which, by 
strongly outcoupling the cavity to a detector (e.g. the 
CBJJ of Fig. [I} the decay rate of the cavity is engineered 
to be much larger than the energy relaxation rate of each 
qubit. 

By continuously monitoring the cavity in a period 
longer than the cavity life time but much shorter than 
the qubit relaxation time, we can know the state of the 
two qubits. If a photon is emitted out the cavity, the sys- 
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FIG. 9: The Q value as a function of time. The solid line is 
for an initial excitation in one of the qubits, and the dashed 
line for an initial excitation in the cavity. 



tern collapses to the ground state |0) = | ft) <£> |0) p . Then 
we start from the beginning by re-exciting one qubit and 
repeating the monitoring. If no photon is detected, the 
system is in the eigenstate | A3) , which means a two-qubit 
entangled state is prepared. 

Experimentally, the monitoring can be done by biasing 
the CBJJ to be in resonance with the resonator, |wio — 
u) r \ <C k, where 



1 ci 3 

2C S3 Ld 



(61) 



is the coupling strength between the CBJJ and the res- 
onator, and the coupling has the form (see Eq. (|A19I) ) 



H 



R-C 



(62) 



With an asymmetric design of the stripline resonator 
[35l | . such that the line has only one coupling capacitor 
C m 3 at one end, the total decay rate of the cavity is 
approximately equal to n. 

In order to make sure that the photon absorbed by the 
CBJJ will never go back to the cavity, the macroscopic 
quantum tunneling (MQT) rate Ti for the excited state 
|1) of the CBJJ should be much larger than k [36[. Thus 
the barrier height Uo of the CBJJ should be close to 
1.5wio, as shown in Fig. [TWa). When a CBJJ is excited 
by absorbing a photon, it immediately switches to the 
dissipative branch and creates a voltage pulse which can 
be measured easily either directly or by the use of an 
additional dc Squid [37l ]. 

The MQT rate for the ground state |0) in this case can 
be calculated as 13811 



r °-2^ 1 



'8647rC/ 



■ exp 



36<7 
5w B 



2 x 10~ Vo, (63) 



by assuming the plasma frequency lo p ~ luiq. The rate 
Ti is about 500 times of r . 

To achieve the desired effect with good enough effi- 
ciency, the cavity decay rate n should be in the range of 

r < k < r x . 



(b) 




FIG. 10: Schematic potential energy diagram of the CBJJ 
(phase regime). 



This method is also suitable for the short resonator 
case. As derived in Appendix [Bl the Hamiltonian for the 
qubits-CBJJ system reads 



tt \ " hjj wio £12 

n = — > er„- (7, H <7i (Tn 



J=l,2 



3=1,2 



If the coupling capacitance C m3 3> C m i^2, the direct cou- 
pling £12 between the two qubits is negligible, and the 
CBJJ itself acts as a lossy cavity. By biasing the CBJJ 
on resonance with the qubits, the Hamiltonian projected 
onto the states |1) = | |t> ® |0), |2) = | Tl) ® |0) and 
|3) = I TT}® |1) is then 



H 









hi 








hi 




-«72 






(64) 



which has eigenvalues and eigenvectors similar to those 
in Eqs. ® and ([57]) . 

In this case, the cavity decay rate is given by the MQT 
rate of the upper level Ti . It has to be much larger than 
the qubit decay rates, but not necessarily to be as large 
as that in the long resonator case. So we can bias the 
CBJJ barrier Uq to be higher than 1.5wio. As illustrated 
in Fig. [TOT b'). when Uo ~ 1.84wio, Ti is approximately 
equal to O.Olwio, which can be of the same order as the 
couplings 7^. 

To study the conditional time evolution of the system 
state, we introduce a non-Hcrmitian Hamiltonian [39( 

-ffcond = H -iT 2J a t a I _ i ^i <7 3 (T 3- ( 65 ) 
3=1,2 

Here we assume that both qubits are equally coupled to 
the CBJJ 71 = 72 = 7, and have the same relaxation 
rate T < Ti, 7. 

The eigenvalues and corresponding eigenvectors of 
-ffcond 9-re 

Ab^-^^rii^-r^, A^ = -ir, (66) 



13 



and 



ia; i2 > 



Ti ± i v / 2~f 2 - T'f - 



+ |2)) + ^|3), 



1 



|A£) = -^=(11)- 



(67) 



For an initial state, say, 

n(|Ai) - |A' 2 }) 



|^(0)) 



1 



- ^ (68) 



the (unnormalizcd) state vector at later time t is given 
by 



\m) 



V^ 2 - sr? 

1 

72 



e-* A i*|Ai) - e- a =*|A 2 ) 



(69) 



as shown in Fig. QT] On a time scale T^ 1 <C £ <C 

the conditional probability of preparing the maximally 

entangled two-qubit state approaches unity. 

In a real experiment the fidelity of the states prepared 
by this procedure will not be exactly 1, due to the de- 
tector inefficiency and to the spontaneous decay of the 
qubits during the time required to reach the asymptotic 
regime (less than 100 ns, Fig. [TTj) . Detector efficiencies 
(visibilities) as large as 0.89 have been obtained recently 
for phase qubits [40[ . The contribution of both these pro- 
cesses to the fidelity of the final state in the asymptotic 
regime can be calculated 33j. With parameters for su- 
perconducting qubits as given in Fig. [TTJ we estimate 
that the fidelity of preparing a Bell state by this proce- 
dure will be as high as 90%. 

We conclude this part of the paper by pointing out that 
this method of producing entangled states is robust with 
respect to spurious excitations in the resonator [4l| and 
that can be readily applied to existing experimental ar- 
chitectures such as phase qubit - resonator - phase qubit 
42] without significant changes in the sample design. 




40 60 

time [ns] 

FIG. 11: The conditional occupation probabilities for the 
eigenstates |Ai 2) (the dashed line) and | A3) (the solid line 
from iJcond and the dotted line from Eq. (|70(l ). We have set 
Eji, 2 = low = 10 GHz, 71 = 72 = 100 MHz, 17 = 50 MHz 
and T = 0.5 MHz. 

Thus the probability for finding the two qubits in the 
maximally entangled state (| T4) — I it))/V2, conditioned 
on that no switching event has been detected, is obtained 
from Eq.® as 



VII. CONCLUSION 

We have analyzed the entangling properties of a sys- 
tem consisting of two superconducting qubits coupled to 
electromagnetic fields, both in the quantum and classi- 
cal cases, and in a variety of experimental situations. 
Efficient decoupling of the qubits, two-qubit entangle- 
ment and high-fidelity quantum gates can be obtained 
in this model. We propose also a quantum nondemoli- 
tion scheme for creating qubit-qubit maximally entangled 
states by monitoring the state of the resonator. Our cal- 
culations are not dependent on the physical realization 
of the qubit and the coupling, therefore most of our re- 
sults can be adapted to various qubit species and coupling 
schemes. 
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V^ 2 " 2T'f 



V2 



(70) 



This expression approximately fits the numerical solution 
calculated by the conditional Hamiltonian in Eq. (|6"S")) . 



APPENDIX A: DERIVATION OF THE 
HAMILTONIAN FOR A LONG RESONATOR 

In this appendix, we derive the Hamiltonian for the two 
SCB charge qubits and the CBJJ capacitively coupled to 
a ID stripline resonator. We assume the length of the 
resonator is L, then take c and I as the capacitance and 
the inductance per unit length. 
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We start with the classical Lagrangian of the resonator. 
At a point < x < L, the local charge density q r (x,t) 
and phase ip r (x,t) satisfy the following relations 



d d(p r (x,t) _ g r (x,t) 
dt c ' 



.(x,t) d d 2 cp r (x,t) 



dt 



Wx 2 



where i9 = h/2e. The Lagrangian of the resonator can 
be written as 



S 2 r 1 



L dx {* 


dip r {x, t) 


2 1 


dip r {x, t) 






I 


dx 



Because of the boundary conditions 



dip r (x, t) 



d(f r (x,t) 



x=0 



Ox 



= 0, 



(Al) 
(A2) 



x—L 



(p r (x,t) has the form <p r (t) cos(k n x), with k n = rnr/L. 
By adopting simplified notations tp r (t) — * <p T an d 
d(p r (t)/dt — > (p r , the Lagrangian can be rewritten as 



dx 



cip 2 cos 2 (k n x) ipjk 2 sm 2 (k n x) 



2/ 



(A3) 



The Lagrangians of the SCBs and the CBJJ are as 



follows. 

^SCB 



j=l,2 3 = 1,2 

+ E Ejjcosifj, (A4) 



= 1,2 



c 

£cbjj = —($ (p3) 2 + Ej cos <p 3 + I b d <P3, (A5) 



with the gate voltages V g j = Vd C j + V ac j . 

By considering that the SCBs are close to the ends of 
the resonator x\ — ► 0, x 2 — > L, the CBJJ coupled to the 
end of the resonator x 3 = L, and taking the mode of the 
resonator n — 2, the Lagrangian for the couplings reads 



E 

3=1,2,3 

EC mi 



■[&Q<p r (xj) + 'da'Pj] 



{$0<Pr + VofioY 



= 1,2,3 

The total Lagrangian is then 

£ = £r + £scb + £cbjj 



C 



coup ■ 



(A6) 



(A7) 



The conjugate charges arc calculated by the Legcn- 
dre transformation Qj = —2erij = d£/d(i)oipj) (j = 
r, 1, 2, 3), and the total Hamiltonian 



J 



H = E Qi^o<Pi - £ 

j=r,l,2,3 

(2e) 2 /C S1 C S2 C S3 \ 



2 V 

\2 



n 



2ir' 



J r 2(2e) 2 V Li 



m 2 , (2e 



/ 1 , fm iCt,2Ct,z\ 



n 



(ni - n s i) 2 - Eji cos (px 



, (2e) 2 ^ , C 2 2 C S1 C S3 ^ , , 2 _ , (2e) 2 / C^ 3 C S1 C S2 ^ 2 _ ft 
+ 1 + tt ( n2 " n «a) " ^2 cosv? 2 + — — 1 + mJ ni - Ej cos<^ 3 - Trhvs 



2C S2 V 

(2e) 2 C m iCi; 2 Cs;3 



2C S3 V~ " J 



n 



(2e) 2 C m2 CsiCs3 (2e) 2 C m 3Cx;iCi; 2 
ra r (m - n s i) n r (n 2 - n g2 j n r n 3 



n 



(2e) 2 C m iC m 3Cs 2 , , (2e) 2 C m2 C m 3Csi . (2e) 2 C m iC m2 Cs3 . .. . 

-(m - n g i)n 3 H (n 2 - n g2 )n 3 H (m - n gl )(n 2 - n g2 ), 

(A8) 



n 



n 



where n g] - = — C g j V g j/2e, and 

n :l ~ c 

with 



OS4 I OeiOe 2 OE3 — C/ ml OS20s3 

/~i2 r~i f i /~i2 f i f i 

"°m2°SlL-E3 — O to3 OsiOe 2 , 



Cm — C\ + C g i + C rn i, Cy,2 — C 2 + C g2 + C m2 , 
CE3 = C + C m3l Cs4 = C m i + C m 2 + C m3 . 



By assuming that C m j -C Cnj {j = 1, 2, 3), and defin- 
ing 

_ La Li 

C*Er = -r- + Cs4 and L Er = 77-77 

2 27T Z 

as the effective capacitance and effective inductance of 
the resonator, the Hamiltonian in Eq. (|A8[) can be rewrit- 
ten as 



H 



(2e) 2 a fi 2 
2C Sr r+ 2(2e) 2 L Er 
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£ 

.3=1,2 

(2e) 2 
(2e) 2 



2C 



K" -™j«) 2 -Ej 3 COS Lp j 



75/ COS ¥> 3 - —4^3 

2e 



C* 



771-3 



Cs3 



7 = 1-2 



(2e) 2 C m iC m 2, x 
~p; ^ n — (™i - %i)("2 - %2j 



(2e) 2 C 



m3 



CsrCs3 .' t -7* L-'Sj 
.7 = 1,2 



(A9) 



For a relatively long resonator Lc ^> Cs4 , the direct SCB- 
SCB and SCB-CBJJ couplings, described by the last two 
lines in Eq. (|A9j) . arc negligible. To obtain the quan- 
tum Hamiltonian, we replace the variables rij , (pj by the 
operators hj, <fij which obey the commutation relation 

[&,»*] = Wife, (j,k = r, 1,2, 3). (A10) 

The quantized Hamiltonian of the resonator is then 



(2e) 2 a 
-n 



2C Er r ^2(2e) 2 L & 



So) r [ a T a + — 



(All) 



where (a) is photon creation (annihilation) operator, 
and the resonance frequency 



U) r = l/ \J LzrC-£r ~ 27T / LV fc 

The dimensionlcss charge operator 



i 

Ye 



(a f -a). 



(A12) 



(A13) 



By projecting the Hamiltonian of SCBs onto the charge 
states \rij), we can get 



{hj-n g iY - Ejj cos ip j 
{2ef 



hscb — 


■(2e) 2 
2Cjy 


= £< 




3=1,2 






2 





2C 



(A14) 



We also assume that the two SCBs are in charge regime, 
Eci,2 = 2e 2 /C S i,2 > Ejx t 2, so that we truncate (|A14|) 
to the two lowest charge states of each SCB, and obtain 
the Hamiltonian 



H 



SCB 



£ 

7=1,2 



E 



■Cj n gj 



2 J 2 



(A15) 



form in Eq. ([2]), by transforming the Hamiltonian (|A15|) 
into the uncoupled eigenbasis {| | j), | ||),| IT), I II)}, 
we arrive at 



Hscb = 



£ 



E C jWj(t) cos(ui d t + 4>j)a x 6 - L 



J 3 n- 



(A16) 

where Wj(t) = -C g3 V flwj (t)/2e. 

For the CBJJ, we consider E c = 2e 2 /C E3 < 75j. It is 
better to discuss it in the "position" space, where <p 3 acts 
as a position operator. In the bottom of one of its local 
minima (see Fig. I10[) . the tilted cosine potential is ap- 
proximated by a harmonic potential. So the Hamiltonian 
is approximately 



ifcBJJ — 



2Cyj3 

fkj 10 . 



E j cos ip 3 - —I b ip 3 
2e 



|0)(0|) 



(A17) 



Here |0) and |1) indicate the ground and the first excited 
states of the CBJJ, not charge states anymore. fi^io is 
the energy difference between the two states. The dimcn- 
sionless charge is analogous to the momentum 



n 3 = 



( CT 3 -°3 ) = 7T 



2(2e) 2 ~ J/ 2e V 2 



/^loCsa^^ (Alg) 



|1)(0| and the lowering 



with the raising operator erg 
operator er^ = |0)(1|. 

With these charge operators, we obtain the quan- 
tized Hamiltonians for resonator-CBJJ and resonator- 
SCB couplings from the fourth line in Eq. (|A9[) 

Hr_ c = k{<4 - <? 3 ){a 1 - a), (A19) 
Hrs = i ^2 9] [(J* + 2Wj(t) cos(u d t + fa)] (a f - a), 



7=1,2 



where the coupling strengths 



(A20) 



foC m3 \2u)\QU) r 



fJj 



c 



£7 



% < A21 > 



2 V C S3 Lc' 
The total Hamiltonian is 

H = TTTr + i?scB + flcBjj + Ar-s + fik-c- (A22) 

For a very short resonator Cy> — * Cyj4 and n r — > 0, 
the kinetic energy terms in Eq. (|A9[) become the same as 
terms in Eq. (|B3[) . furthermore, ftuv is much larger than 
the other energies in the system, therefore we can assume 
that the two SCBs and the CBJJ are directly coupled to 
each other. 



APPENDIX B: DERIVATION OF THE 
HAMILTONIAN FOR A SHORT RESONATOR 



Since the DC gate voltages are biased at the charge 
co-dcgcneracy point (TT|) and the AC voltages have the 



When the stripline resonator is very short, Lc -C Cyj4, 
the center conductor can be considered as an island. The 
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SCBs and the CBJJ are capacitivcly coupled to this is- 
land. Instead of deriving the Hamiltonian from the clas- 
sical Lagrangian of this circuit, here we use a relatively 
simpler method. 

The system has four nodes. As shown in Fig. [IJ the 
total charge on node j (j = 1,2,3) is indicated by Qj = 
—2erij (2 for a Cooper pair). The center island acts as 
the fourth node, and the total charge on it is assumed to 
be Qi = — 2eri4. Since the total charge on a node is equal 
to the sum of the charges on all capacitors connected to 
the node, by denoting the electrostatic potential of node 
j as Vj, we write the charge equations for all the nodes 
in a matrix form, as 



Qi 


~ CglVgl ' 




' Cei 








Cml 




\Vi 1 


Qi 


- G g2 Vg2 












C m 2 




v 2 




Q 3 












C m 3 




v 3 








- Cml 


C m 2 


Cm 3 


Cs4 . 




.v 4 _ 



(Bl) 

The 4x4 matrix in Eq. (|B1|) is called the capacitance 
matrix C. The total electrostatic (kinetic) energy of the 
system can be calculated by using the matrix 



T 



(B2) 



where Q denotes the column vector of charges on the left- 
hand side of Eq. (|B1[) . By making the same assumption 
as in Appendix IA1 C m \.2,3 <C Cei, 2, 3, and take Q4 = 0, 
the result of Eq. (JB2J is ' 



T 



E C \{ni - n g i) 2 + E C2 (n 2 - 
+[Eis(ni - n g i) + E 23 (n 2 
+E 12 (ni - n g i)(n 2 - n g2 ), 



n g2 ) + 
- n g2 )]n 3 



Ecnj 



(B3) 



with 



E, 



ci 



E, 



C2 



E f 



* 2e 2 /C sl 
a 2e 2 /C E2 
2e 2 /C S3 , 



E t 



E 



23 



En 



(2e) 2 C m iC m3 

S~1 SI S~1 ' 

OsiL/S3L'S4 

{2e) 2 C m2 C m3 

Cs2C , S3Cs4 

{2e) 2 C m iC m2 

C-£lC-£ 2 C-£4 



(B4) 



The inductive (potential) energy of the system can be 
expressed as 



U = -Ejt cos ipx - E J2 cos ip 2 - Ej cos ip 3 - —I b ip 3 , 

2e 

(B5) 

where ipj (j = 1,2, 3) is the gauge- invariant phase differ- 
ence across each Josephson junction. 

The total classical Hamiltonian is H = T + U. To 
derive the quantum Hamiltonian, we again replace the 
variables ipj and nj by operators if>j and fij. 

The quantized Hamiltonians of the SCBs and the CBJJ 
have the same form as those in Eqs. (|A16[) and (IA17I) . 
With the charge operators derived in Appendix [AJ the 



Hamiltonian for direct coupling between the two SCBs is 
given by 

-ffs_s = E 12 (hi - n g i)(h 2 - n g2 ) 

H\2 x x , E\2 , , , . , 1 \ x 

= —&fv 2 + —W2[t) cos(uj d t + <\> 2 )<jI 
+-^t«i(t)cos(w (J t + 0iK, (B6) 

and the Hamiltonian for the couplings between the two 
SCBs and the CBJJ is derived from the second line in 
Eq. 



s-c 



where 



3=1,2 



lj = 2e 



a j a 3 + w j(t) cos(uj d t + 4>j)al 



C 



Cr, 



(B7) 
(B8) 



Cs4 \ Cjy ) V 2Ce3 

The total Hamiltonian of this system is finally 

H = i?scB + -ffcBjj + Hss + Hs-c- (B9) 



APPENDIX C: DERIVATION OF THE 
EFFECTIVE HAMILTONIAN IN THE 
DISPERSIVE REGIME 

We assume that the Rabi frequencies ilj (t) change adi- 
abatically with respect to the qubit (Larmor) frequency, 
so in the following calculations we take dSlj/dt = 0. We 
use the Baker-Hausdorff formula to expand the transfor- 
mation to second order 



H A 



1 f) A 

H+[H,A] + -[[H,A},A}- i — 



E 



1 



E 

i=l,2 ]UJr 

- VLj cos(u!d + (f>j)<Jj + u) r a) a 
3=1,2 

.91.92 ( J_ J_ 



(a+a 2 -+ara 2 + ). (CI) 



By making an adiabatic approximation for the res- 
onator, (a^a) = and (a' + a) = 1, we neglect the ac- 
Stark shift and Lamb shift terms and the Hamiltonian of 
the resonator. The effective Hamiltonian can be written 

MS 



Ha = Hao + Hai, 



where 



H A0 = E 

3=1,2 



^""amiwdt + hW, (C2) 



AjLU r 
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and 



APPENDIX E: QUBIT-TLS SYSTEM 



H A1 = £ 

.7=1,2 

5i32(Ai + A 2 



^-cr| + % cos(w d i + 0j)cr? 



(at a 



2 



(C3) 



2A!A 2 

Now we transform Hai into the interaction picture 
jj = e «/o dt ' H Ao(t')jj Aie -iJo dt'H A0 (t')_ 

The raising and lowering operators become 



]T *"J„(±z)e 



(C5) 



where J„(±z) is the n-th order Besscl function of the first 
kind. Since z = gjClj j (AjU) r ) <C 1, Jo(±z) is dominant. 
Therefore ~ , and the Hamiltonian 



H 



j=l,2 



/.7 



; cr| + % cos(o> d i + <^-)cr? 



2AiAz W o 2 + cr 1 cr+). (C6) 



A lot of experimental progress has been made recently 
on phase qubits following the realization that the dielec- 
tric insulator forming the Josephson junction contains 
two- level system (TLS) defects [3, [32[. These defects 
have been shown to have decoherence times comparable 
occupation prob. of |2) 




0.03 0.04 0.05 0.06 0.07 0.08 0.09 

£2/ Ei 



FIG. 12: (Color online) Leakage for a 7r-pulse. 



APPENDIX D: LEAKAGE OUT OF THE 
COMPUTATIONAL SUBSPACE 



Here we only consider a single charge qubit with rect- 
angular driving pulse f2 cos(wrft). By considering the 
three lowest charge states |n) = {|0), |1), |2)}, and dc 
biasing to the degeneracy point, the Hamiltonian has the 
form 



H 



E c (n - n g ) 2 




— Ej cos ip 

-Ej/2 
-Ej/2 -20cos(w d *) 
-Ej/2 



2Er 





-Ej/2 
— 4f2 cos(ojdt) 



The leakage can be defined as the occupation probabil- 
ity of state |2) after certain pulse duration, averaged over- 



all possible initial states j^i, 



10) 



sin § 1 1) 



In Fig. [121 the leakage for a 7r-pulsc is plotted by nu- 
merically solving the Schrodinger equation with H . The 
driving frequency is set to u>d = Ej. 



to that of the qubit, thus they can be addressed coher- 
ently (e.g. by tuning the qubit on- and off- resonance 
with them). 

The form of the interaction Hamiltonian between the 
qubit and the TLS is of the type a x a x in the case of phase 
qubits [32|, and this coupling becomes important when 
A = \w\ — UJ2\ < oj xx . Here we adopt the same notations 
as in Sec. Mil u> xx denotes the coupling strength between 
the qubit and the TLS, u>\ and w 2 are Larmor frequencies 
of the qubit and the TLS, respectively. By assuming 
that for a single qubit there is only one such TLS near 
it, and the TLS is weakly coupled to the driving field 
of the qubit, we may use the Hamiltonian in Eq. (|24|) 
to describe this qubit-TLS system, and therefore use the 
switchable scheme developed in Sec. MI Al to perform 
quantum gates with the qubit and the TLS. 
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